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ABSTRACT 
The Hypervelocity Star survey presents the currently largest sample of radial velocity measurements of halo 
stars out to 80 kpc. We apply spherical Jeans modeling to these data in order to derive the mass profile of the 
Galaxy. We restrict the analysis to distances larger than 25 kpc from the Galactic center, where the density 
profile of halo stars is well approximated by a single power law with logarithmic slope between -3.5 and 
-4.5. With this restriction, we also avoid the complication of modeling a flattened Galactic disk. In the range 
25 < r < 80 kpc, the radial velocity dispersion declines remarkably little; a robust measure of its logarithmic 
slope is between -0.05 and -O.L The circular velocity profile also declines remarkably little with radius. 
The allowed range of Ve(80kpc) lies between 175 and 231 km s~ , with the most likely value 193 km s~ . 
Compared with the value at the solar location, the Galactic circular velocity declines by less than 20% over an 
order of magnitude in radius. Such a flat profile requires a massive and extended dark matter halo. The mass 
enclosed within 80 kpc is 6.9;!;[ " x 10" M0. Our sample of radial velocities is large enough that the biggest 
uncertainty in the mass is not statistical but systematic, dominated by the density slope and anisotropy of the 
tracer population. Further progress requires modeling observed datasets within realistic simulations of galaxy 
formation. 
Subject headings: Galaxy: formation — Galaxy: halo — Galaxy: kinematics and dynamics 



1. INTRODUCTION 

Measuring the mass of the Galaxy is an astonishingly diffi- 
cult task. Convenient tracers of disk rotation - stars an d gas 
clouds - extend only to 20 kpc fe.g. ISofue et al.ll2009h . At 
larger radii, statistical analysis of radial velocities must be 
used. Traditional tracers at distances up to 100 k pc include 
globular cluster s and dwa rf satellit e galaxies (e.g..lKochanekl 
Il996; Wilkinso n & Eyansil999: Sakamoto etal.ll2003h . Re- 
cently, Battagli a et al.l (l2005h derived the radial velocity dis- 
persion profile to 120 kpc using a combined sample of glob- 
ular clusters, satellite galaxies, and halo red giant stars. They 
found the dispersion falling from ~ 1 20 km s~' t o^50 km s~' 
between 20 and 120 kpc. In contrast. fXue et al] (l2008h assem- 
bled a large sample of blue horizontal-branch (BHB) stars 
from the Sloan Digital Sky Survey Data Release 6 (SDSS 
DR6) and found a much flatter profile bet ween 20 and 60 kpc . 
Here we use a new spectroscopic survey jBrown et al.ll20ld) 
aimed at finding hypervelocity stars (HVS) to set the most 
p recise constraint on the Galactic mass within 80 kpc. 

iBrown et al.l (120101) present a sample of 910 late B-stars and 
early A-stars in the Galactic halo. Their luminosity, and there- 
fore distance, depends on whether these stars are BHB stars or 
main-sequence blue stragglers with similar effective temper- 
ature and surface gravity. The ambiguous nature of the stars 
is especially problematic at redder colors, u-g > 0.6, where 
the lumino sity can differ by a factor of 5, as demonstrated 
in Fig. 4 in lBrown et al.l (I2010I) . This bimodal distribution of 
distance does not have a well-defin ed average value and there- 
fore requires statistical sampling. iBrown et al.l (|2010) create 
100 Monte Carlo realizations of each star to sample the color 
and metallicity distributions and to derive the distributions of 
luminosity and distance. We use this full Monte Carlo cata- 
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log of distances in our analysis, while retaining the observed 
values of radial velocity and its uncertainty. 

2. METHOD 



We use the spherical Jeans equation (iBinney & Tremaina 
I2OO8) to solve for the circular velocity Vc given the radial ve- 
locity dispersion a,- of tracer particles: 
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where M{r) is the enclosed total mass within radius r, ptr is the 
density of tracer population, and /3 is the anisotropy parame- 
ter. This equation assumes that the mass distribution M{r) is 
static and spherically symmetric. 

In our case the tracers are halo stars, which follow a steep 
density profile with negative logarithmic slope (see below) 
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The velocity dispersion slope and the anisotropy parameter 
are typically less than unity, and therefore, 7ti dominates the 
last factor in equation ([U. The circular velocity profile de- 
pends mainly on 0-^ and 711. 

For ease of comparison with mass profiles expected from 
cosmological simulations, we fit a power law relation to the 
radial velocity dispersion, normaUzed at radius ro: 



ar(r) = ao (r/ro) 



(3) 



We apply a maximum like lihood (ML) method , similar to that 
described in Appendix of iGnedin et al.l ilOOT), to maximize 
the probability of the model fit given the observations: 
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TABLE 1 

Characteristic Parameters of the Mass Profile 
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Fig. 1 . — Density profile of tracer populations: BHB stars from HVS sur- 
vey (solid circles), BHB stai's from SDSS DR6 (open circles), and RR Lyrae 
from SDSS Stripe 82 (solid squares). The density is calculated directly froin 
observed counts, without correcting for sky incompleteness. The true density 
is a factor ~ 5 higher Errors are from Poisson statistics. Lines show best fits 
for the logarithinic slope, from 3.5 (dashed) to 4 (solid) to 4.5 (dotted). 



We assume that each radial velocity Vrj is drawn from a Gaus- 
sian distribution with zero mean and with combined variance 
of the intrinsic dispersion at its location, o"r(r,), and the mea- 
surement uncertainty ct,,_,-. The ML method uses all available 
observational information without binning the data. 

The value of the velocity dispersion is sensitive to the pres- 
ence of outliers in the sample. In order to remove them, we 
impose a cut |v,..,| < Vesc, with the escape velocity calculated 
self-consistently from the derived mass profile. In the "most- 
likely" case this process removes only 6 out of 558 stars lo- 
cated at r > 25 kpc. 

The escape velocity at a given radius r measures the total 
gravitational potential, a sum of the two terms interior and 
exterior to r: 



y2 =. 

' esc — 



-2$(r) : 
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We evaluate the second term by adopting a power-law profile 
for the total matter density pir) oc r"^"", out to a maximum 
radius rout- In the outer Galactic halo we expect < a < 1. 
Then the escape velocity is 
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Ignoring the boundary condition at rout, we obtain 
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IWaikins et al.l (120101) show that most halo tracers reside in the 
radial range where a w 0.5, and thus yesc('") ~ 2y.(r). As 
the resulting circular velocity profile is close to flat, we adopt 
this simple relation in our analysis. The ML fit converges to 
the same value regardless of the initial guess for ygsc- In the 
inner regions of the Galaxy, lower a leads to higher escape 
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velocity, consistent with the estimate from the RAVE survey 
ofJ4sc ~ 500-600 km s~' at the solar location JSmith et"al] 
l2007h . 

The density profil e of RR Lyrae halo stars shows a brea k 
around 25-30 kpc dWatkins et alj 120091: ISesar et all lIoToh . 
with shallow inner slope (7;^ « 3) and steep outer slope (7ti « 
4.5). The stellar halo is also less flattened in the outer regions. 
In order to obtain a robust measure of the mass profile, we 
restrict our analysis only to radii r > 25 kpc. Figure [T] shows 
that at these radii the density profiles of both BHB stars and 
RR Lyrae are consistent with our fiducial value 7tr « 4. 

We consider the following parameter ranges: 3.5 < 7tr < 
4.5 and < /? < 0.5. For each parameter set, we use the 
full Monte Carlo catalog and self-consistently remove outliers 
using the value of the escape velocity resulting from the ve- 
locity dispersion profile. This procedure gives simultaneous 
estimates for both ar(r) and Vdr). 

3. RESULTS 

The radial velocity dispersion profile declines remarkably 
little with radius. For the power-law fit (eq. |[3l) we obtain 
(To « 1 14 km s~' and 7cr ~ 0. 1 for ro = 40 kpc. This choice of 
ro minimizes the error of (Tq- We do not quote formal errors of 
the fit because the uncertainty is dominated by the unknown 
parameters 7tr and (3, as we discuss below. Table [T] summa- 
rizes the parameters of the best-fit and the most extreme al- 
lo wed models. 

iBrown et al.l (120101) consider a linear fit to the same data and 
obtain ar « 120-0. 3r, where r is in kpc and c, is in km s~'. 
By linearly expanding equation (O near ro, we derive a similar 
e xpression in the ra nge 25 < r < 80 kpc: Ur ~ 120-0. 22r. 

IXue et all (120081) fit an exponential function to their SDSS 
DR6 sample and obtain 0-^ = 111 exp(-r/354kpc) km s~' . A 
linear approximation gives (Jr ~ 111- 0.31 r, which has a 
similar slope and slightly lower normalization than the HVS 
sample. Since our dispersion profile is calculate d from the 
Monte Carlo catalog, we similarly bootstrap the IXue et all 
(l2008i) sample 30 times and obtain a distance distribution, 
from which we calculate the dispersion plotted on Fig.|2] 

The profiles from both samples are less steep than that de- 
rived by £attagliaeLaL(2QQ5J: a,- « 132-0.6r. 

Figure |2] shows the velocity dispersion profiles for all three 
samples. They appear consistent with each other within the 
errors, and the differences in the derived best fit parameters 
can be attributed to small sample sizes. Outside 25 kpc, the 
Battaglia, Xue, and Brown samples contain 80, 741, and 558 
objects, respectively (using the distance distribution from the 
Monte Carlo catalogs). Outside 50 kpc, the three samples 
contain 24, 76, and 163 objects, respectively. Note that our 
sample more than doubles the number of distant stars and thus 
presents the most accurate current measurement of ar- Indi- 
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Fig. 2. — Radial velocity dispersion profile of BHB stars from HVS sur- 
vey (solid circles), from SDSS DR6 (open circles, offset by 1 kpc for 
clarity), and from a combined sample of globular clusters, satellite galax- 
ies, and halo red giants (triangles). Lond-dashed line is a power-law fit, 
<j(r) =111 (r/40)"'"'^ km s"' . Solid line shows the mean streaming radial 
velocity for the HVS sample. Shaded regions illustrate the allowed range 
of circular velocity when the anisotropy parameter is varied from /3 = to 
/3 = 0.5 for a fixed tracer density ju = 4 (bottom left towards top right) and 
when the tracer density is vaiied from 7tr = 3.5 to ju = 4.5 for a fixed /3 = 0.4 
(top left towards bottom right). Middle solid line i s for 7tr =4, = 0.4. Filled 
circles connected by vertical lines show the Watk ins et alj 120101) mass esti- 
mator applied to the HVS sample at 40, 60, and 80 kpc. 



vidual velocity errors and resampling of the distance distribu- 
tion affect the value of (7^ by less than 3%. The streaming ve- 
locity (vr) is small everywhere except at the outermost radii. 
Overall, the data indicate that the velocity dispersion varies 
little with radius out to 80 kpc. 

We can now combine equations dU and ^ to calculate the 
circular velocity 



Vc(r) = ao {r/ro) ~^' (7„ + 27, - 2/3) 



1/2 



(8) 



This estimate is degenerate with respect to the density and 
anisotropy of the tracer population. Lower 7tr and higher /3 
result in the lower mass estimate, and vice versa. At the same 
time, steeper density slope can balance stronger anisotropy. 
Based on our derivation of the tracer density in Figure [T] we 
believe the slope is constrained to be between 3.5 and 4.5, 
with the most likely value of 7ti w 4. The anisotropy pa- 
rameter is not directly known, but we can take the predic- 
tions of cosmological simulations of galaxy formation as a 
guide. The centers of halos are close to being isotropic. In 
the outer parts, the orbits of dark matter particles and satel- 
lite halos become more radially biased with distance, reaching 
P Ki 0.5 out side the peak of the circular v elocity curve (e.g., 
lDiemand"et al. 2007; Nav arro et al.ll2010l) . Thus a full range 
of values < /3 < 0.5 is possible, while the most likely value 
in our radial range is /3 « 0.4. 

Figure |2] illustrates the allowed range of circular velocity. 
We take the best fit {^^ = 4, /3 = 0.4) and consider the full 
range of variation of either parameter while keeping the other 
fixed. The overlap of the two shaded regions gives the most 
probable values of Vc- The uncertainty due to either parameter 
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Fig. 3. — Shaded regions illustrate the allowed range of circular velocity 
when the anisotropy parameter is varied from /9 = to /3 = 0.5 for a fixed 
tracer density 7(r = 4 (bottom left towards top right) and when the tracer den- 
sity is varied from 71, = 3.5 to7i, = 4.5 for a fixed /3 = 0.4 (top left towards 
bottom right), same as in Fig.|2] Thick dashed line shows a representative fit 
for an NFW halo model (see text for details). Independent const raints at 8, 50, 
60, and 100 kpc are fro m M cMillan & Binney 12 010), Wilki nson & EvansI 
fTg99ll , PCue"etarU200g) . and lWatkins eTaP JMOl) , respectively. 



is systematic in nature and comparable in magnitude. The 
value of the circular velocity at 80 kpc is uncertain at least by 
10%, and the value of the enclosed mass at least by 20%. The 
most likely values are ^(80) = 193 km s"', M(80) = 6.9 x 
10" Mq. 

For comparison, w e also use an indepen dent robust mass 
estimator proposed by lWatkins et al .W2010). The mass within 
radius r is given by the following average of radial velocities 
of all objects inside r. 



V>) = (7tr + a-2/3)(v?,(^)"), 



(9) 



where (-2 - a) is again the logarithmic slope of the density 
profile in the range of radii probed by the data. This estimator 
also depends on 7tr and /3, similarly to equation ([8]l. In fact, 
if these two parameters and a are constant with radius, then 
a = 27ct. The estimator is different in that it takes a single 
average of all radial velocities instead of fitting the dispersion 
profile. 

Three vertical lines in Figure |2] show the applications of the 
robust estimator at 40, 60, and 80 kpc. They indicate a slightly 
steeper decline of Vc with radius, but still fully consistent with 
the results derived from the dispersion. The most likely value 
at 80 kpc is V,,^,, (80) « 190 km s"'. 

Figure |3] compares our derived circular velocity with other 
independent observational constraints. Dotted lines show the 
extrapolation of Vc outside the range of our data and should be 
treated with caution. Watkins et al. (2010) estimate the mass 
of the Galaxy at r = 100 kpc using the sample of all known 
satellite galaxies. Depending on the anisotropy and the inclu- 
sion of particular galaxies, the mass can vary between 3.3 and 
13.8 X 10" Mq, corresponding to the circular velocity range 
from 1 19 to 244 km s~'. Our results favor the upper half of 
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this range. 

Within 60 kpc, IXue et alj (12008b estimate K = 170 ± 
15 km s~' by matching observed velocity dispersion to the 
motion of particles in two hy drodynamic simulations o f 
galaxy formation. Within 50 kpc, [Wilkinson & Evani ( 119991) 
obtain M(50) = 5.4:!:3g x 10" Mq by modeling the distribu- 
tion function of observed veloci ties. For c l arity, we do not 
show similar results obtained by Kochanek ('1996'), M(50) = 
4.9!|j X 10" Mq, and Sakamoto et al. (2003), M(50) = 
5.5!!]:{x 10" Mq. 

The circular velocity at the solar circle, Vq, is constrained 
better but it scales with the distan ce to the Galactic cen- 
ter, Rq. iMcMillan & BinnevI (12010 ") analyze the motion of 
masers in star-forming regions throughout the Galaxy and 
find that the best-constrained parameter is the ratio Vq/Rq = 
30 ± 2 km s~' kpc"' . The distance Rq can vary from 6.7 to 8.9 
kpc, depending on model assumptions. We plot this constraint 
on Vb as a diagonal line in Figure[3] This line matches the ex- 
trapolation of our fits perfectly. However, this result should 
not be considered as an improvement on the value of Vq, be- 
cause we do not expect Vdr) to remain a power law with fixed 
slope at such small radii. 

For illustration, we also show an example of a three- 
component model of the Galactic potential, sim ilar to the 
cosmologically-motivated model bv lKlvpin et al.l (2002). We 
assume an axisymmetric exponential disk with the mass 5 x 
10'° Mq and scale length 3 kpc, and a compact bulge with the 
mass 5 x 10^ Mq dBinnev & Trem aine'2008'). We include a 
dark matter halo represented by an NFW profile with the scale 
radius r, = 20 kpc and the mass that, along with the disk and 
the bulge, gives the total virial mass Mvi,- = 1.6 x 10'^ Mq. 
This model lies in the middle of our inferred interval of cir- 
cular velocity and can be taken as a good first approximation 
to the Galactic mass distribution. Note however that the un- 
certainty in the value of the virial mass is at least 20% and 
possibly larger 

In Figure 2] we illustrate further the dependence of the de- 
rived value of the circular velocity at 80 kpc on the parameters 
7tr and /3. For each combination of the parameters, we run a 
ML fit to the velocity dispersion profile with self-consistent 
removal of unbound stars. The plot shows, however, that con- 
tours of Vc (80) lie close to the lines j^ - 2/3 = const, as ex- 
pected from equation ([8]) for a constant 7^. Larger /? gives 
smaller velocity, but according to cosmological simulations it 
should not exceed 0.5 and be near 0.4. However, the estimate 
of Vc would be higher if the tracers had tangential anisotropy, 
/3 <Q. A density slope 7tr « 4 is most likely, as it fits both the 
BHB and RR Lyrae samples. 

4. SUMMARY 

Using maximum-likelihood analysis of a new sample of 
radial velocities of distant halo stars, we infer that their ra- 
dial velocity dispersion profile declines little with distance 
from the Galactic center in the range 25 < r < 80 kpc: 
a(r) = 111 (r/40kpc)"'''^**km s~'. Spherical Jeans modeling 
indicates that the circular velocity profile Vdr) also falls only 
slightly over the same radial range and reaches between 175 
and 231 km s~' at 80 kpc. The corresponding enclosed mass 
M(80) is between 5.7 x 10" Mq and 1.0 x 10'^ Mq. A 
three-component model for the baryon and dark matter mass 
distribution gives the total virial mass of the Galaxy M^t = 
(1 .6 ± 0.3) X 10'2 Mq at the virial radius T^vir = 300 kpc. 
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Fig. 4. — Contours of the circular velocity at 80 kpc as a function of the 
anisoti'opy pai'ameter /3 and logarithmic slope of the density of tracers 7ir. 
The values of Vc (80) at each point are calculated for the ML fit to the velocity 
dispersion profile. 



Our inferred mass of the Galaxy is higher than that ob- 
tained by Battagha et al.1 (l2005l) (Mvir « 0.8 x lO'^ Mq) and 
IXue et all dlOOSl) (Mvir ^^ 1.0 x 10'^ Mq) based on the mod- 
eling of their radial velocity datasets. Our HVS sample con- 
tains more objects at r > 40 kpc than the Battaglia et al. and 
Xue et al. datasets. Thus we have a stronger constraint on 
the shallow slope of the velocity dispersion profile and we 
derive a correspondingly larger mass. Our infe rred mass is 
consistent with the larger scale measurement by iLi & Whita 
(2008) based on the Andromeda-Milky Way timing argument, 
Mvii- ~ 2.4 X 10'^ Mq. The implied dynamical mass-to-light 
ratio of the Galaxy, M^a/Ly ~ 50 in solar units, is also con- 
sistent with galaxy-galax y weak lensing measurements by 
'Mandelbaum e t all (120061) . galaxy kinematics modeling by 
More et al. (20 lW, an d halo abundance matching modeling 
by Mo steretal.T(l2010l) . 

Our sample of radial velocities is large enough that the 
biggest uncertainty in the mass estimate is not statistical but 
systematic. Within the framework of spherical Jeans mod- 
eling, the uncertainty is dominated by the density slope and 
anisotropy of the tracer population. These parameters could 
be better constrained by future all-sky surveys of halo BHB 
stars. Deeper surveys that target more distant stars at r > 100 
kpc would be similarly dominated by uncertainty over the un- 
derlying distribution of the tracers. 

The validity of spherical Jeans modeling is also limited 
by the presence of structure in the distribution of halo stars. 
Galactic stellar halo contains remnants of disrupted satellite 
galaxies, some of which are still detectable as tidal streams. 
Stars at ^ 100 kpc from the Galactic center may not have had 
enough dynamical times to reach dynamical equilibrium, fur- 
ther limiting the application of equilibrium modeling. A first 
step in the direction of circumventing these systematics was 
taken by IXue et aL (.200 8,) . who modeled the motion of par- 
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tides in realistic simulated halos. Extension of such analysis 
to many different halo realizations using large samples of ob- 
served velocities may reduce the uncertainty over the global 
mass distribution in the Galaxy. 
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